Rigidity and intermediate phases in glasses driven by speciation 
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The rigid to floppy transitions and the associated intermediate phase in glasses are studied in 
the case where the local structure is not fully determined from the macroscopic concentration. 
The approach uses size increasing cluster approximations and constraint counting algorithms. It is 
shown that the location and the width of the intermediate phase and the corresponding structural, 
mechanical and energetical properties of the network depend crucially on the way local structures 
are selected at a given concentration. The broadening of the intermediate phase is obtained for 
networks combining a large amount of flexible local structural units and a high rate of medium 
range order. 

PACS numbers: 61.43.Fs-61.20.-x 



Concepts from mean-field rigidity in networks found 
their origin from Lagrangian constraint counting in me- 
chanics [l| and have been applied with great success in 
glass science for several decades . Bonds in a glass net- 
work can indeed be considered as constraints arising from 
interatomic stretching and bending forces. The connec- 
tivity or cross-link density (best quantified by the net- 
work mean coordination number f) plays therefore a key 
role. In highly cross-linked networks where r is large, 
there are more constraints than degrees of freedom per 
atom on average and the structure is stressed rigid (hy- 
perstatic or overconstrained). At low connectivity, one 
has a flexible (hypostatic or underconstrained) structure 
that contains more degrees of freedom than constraints. 
Thorpe |3| analyzed the vibrational behaviour of such 
kind of networks and identified a mean-field (MF) floppy 
to rigid transition when the mean coordination number 
equals f — fc — 2.38, a result that agrees with global 
(Maxwell) constraint counting as enunciated by Phillips 
|j| from the enumeration ric of bond stretching and bond 
bending forces. 

The underlying nature of this peculiar transition has 
been deeply reinvestigated recently because two transi- 
tions at ^iid fc(2) have been found Q experimen- 
tally in a variety of glasses. These define an intervening 
region (or intermediate phase, (IP)) between the floppy 
and the stressed rigid phase. In the IP, glasses display 
some remarkable properties such as absence of ageing |y| 
or stress 0, selection of isostatically rigid local struc- 
tures or weak birefringence The two boundaries 
have been characterized from numerical calculations 
and cluster analysis on self-organized networks and 
identified as being a rigidity transition at low r and a 
stress transition at high f. In the mean-field approach 
or in random networks where self-organization does not 
take place, both transitions coalesce into a single one. 
Moreover, links between IP and protein folding pTj |. 
high-temperature superconductors |l2l | or computational 



phase transitions [l3j have been stressed that go much 
beyond simple analogies. The understanding of the IP 
is therefore of general interest. It has become clear that 
stress avoidance in the network is responsable for the 
and the location of the interme- 
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diate phase, an idea that has gained some strength from 
energetical adaptation in a simple random bond model 
|Q| for the rigidity transition or suppressed nucleation 
of rigidity duringa fluid-solid transition . Mousseau 
and co-workers [ig have also shown recently that self- 
organization with equilibration on diluted triangular lat- 
tices would lead to an intermediate phase. 

However, the recent discovery of an IP in more com- 
plex glass systems such as silicates ^3 raises a new chal- 
lenging issue. On the experimental side, most of the re- 
sults have been obtained up to now on simple network 
glasses (e.g. the archetypal GexSei-x), where r can be 
directly related to the concentration [f = 2 + 2x) of the 
species involved 0|. This happens to be not the case 
any more in multicomponent systems such as (even sim- 
ple) binary glasses (a network former, e.g. Si02 and 
a modifier, e.g. Li20) where a non-trivial speciation 
can appear depending on the nature of the cation or 
the atoms involved. This contributes to f in a non- 
linear fashion and application of constraint counting al- 
gorithms becomes more difficult. We describe in this Let- 
ter how speciation affects onset of rigidity and the inter- 
mediate phase. A simple model to highlight the effect 
of the speciation is solved, and combined with cluster- 
constraint calculations applied to silicate or thiosilicates 
of the form (1 - x)SiX2 - XM2X with {X = O, S, Se) 
and {M = Li, Na, K) which are known 19] to display 
a mean-fleld rigidity transition at a; = 0.20. The pres- 
ence of an intermediate phase is demonstrated and its 
structural, mechanical and energetical properties char- 
acterized. It appears that the selection of more flexible 
local structural units with addition of a modifler leads 
to a broadening of the IP, independently of the degree 
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Fig. m with respective probabilities p^l\ pip and 
The behaviour of the pf^^s with modifier concentration 
X can be determined from the normaUzation condition 
pP + Pp + P2^^ — li the charge conservation law 
R = pP + 2pP , and finally the definition of the equi- 
librium constant |2^] of the chemical reaction given 
by: K f, = p^ p^ / p'P p^'-p ■ For instance, in lithium sil- 
icates |25j Ke is of the order of 0.3 at x = 0.17. With 
these equations, the speciation is fully determined with 
respect to x and given by: 

i?(2 - R) 



.(1) 



FIG. 1: An example of a network made of 23 tetrahedra 
with 4 Q^, 4 Q'^ and 15 units, corresponding to the mod- 
ifier concentration x — 0.207. The fraction of edge-sharing 
tetrahedra is rj = 0.32. 



of medium range order. The latter contributes however 
to the increase of Af as well. On the other hand, the 
possibility of the system to adapt the speciation in order 
to lower the constraint free energy in the stressed rigid 
phase leads to a situation that ressembles very much to 
sodium silicates. Taken together, these results provide 
new benchmarks to study IP's in multicomponent oxide 
or chalcogenide glass systems such as fast ionic amor- 
phous superconductors where the speciation and hence- 
forth the elastic nature of the network crucially determine 
physical and electric transport properties. 

In binary sodium silicates (X=0, M=Na), speciation 
depends weakly ^2d\ on the nature of the modifier cation 
M which creates almost only Q"=^ units at low x so that 
the probability of finding the latter is i? = 2a;/(l — x). 
Here, the superscript "n = 3" denotes the number of 
bridging oxygens on a S'i04_„/2M4_„ tetrahedron (a Q" 
unit) that connects to the rest of the network (Fig. 
This means that the chemical reaction l2ll: 



2Q3 ^ q4 ^ q2 



(1) 



involving the species Q" is desequilibrated on the left 
side. A radically different situation is encountered in 
systems with modifier cations of smaller sizes (M=Li) or 
in thiosilicates (X=S) and seleniosilicates (X=Se) [2j|. 
Certain of these glass networks can be indeed made out 
of and species only or, at least, of a mixture of 
all (3"'s. Noteworthy is the fact that Maxwell constraint 
counting does not distinguish between the aforemen- 
tioned systems, although the location of the rigidity and 
stress transitions should be obviously changed. 

Size increasing cluster approximations (SICA) can be 
used |0| to infer the effect of speciation on the loca- 
tion and the width of the intermediate phase. We con- 
sider a network of N tetrahedra Q"*, and (see 
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out of which is obtained = (-R ^ P3 )/2 and pl ' = 
1 — ^3^' ^Pp ■ Starting from this short-range order distri- 
bution pP (the basic SICA units at the initial step I = 1 
which will serve as building blocks), one constructs the 
12 possible structural arrangements of two basic units 

{I = 2), i.e. - Q\ Q4 _ g3^ Qi _ Q2^ Q3 _ Q3^ 

etc. Three energy gains, Egtress, Ei^o and Efiex with 
corresponding Boltzmann factors e; — exp[— Si/ZcsT], 
are defined following the mechanical nature of the cre- 
ated cluster (stressed rigid, isostatically rigid and flexi- 
ble). The probabilities of the created clusters {I = 2) are 
then given by p^j^j oc WkjP^i^''p'j^''ei {i = stress, iso, flex) 
where Wkj is a statistical factor taking into account the 
number of equivalent ways to connect two (^ = 1) units 
together. For instance, there are W44 = 72 different ways 
to connect two Q"* tetrahedra by edges. One should also 
note that there are only two stressed rigid clusters cre- 
ated (corner and edge-sharing Q"* — Q'* connections) and 
one isostatically rigid cluster (rif, = 3.0, a corner-sharin; 
Q** — connection). Maxwell constraint counting 
is then applied on the set of {I = 2) clusters that leads 
to the number of floppy modes of the network given by: 
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where nc{kj) and N^j are respectively the number of me- 
chanical constraints and the number of atoms of the clus- 
ter with probability Pj^- ■ Once this is set and start- 
ing from a flexible (floppy) network where stressed rigid 
dendritic — connections are absent and decreas- 
ing a;, one can investigate at which concentration 2:^(1) 
the network will have a vanishing of Z^^) (rigidity tran- 
sition) and at which concentration Xc(2) the network will 
not be able to avoid stressed rigid dendritic clusters (i.e. 
corner-sharing Q^) any more (stress transition). These 
corner-sharing Q'^'s contribute to percolation of stressed 
rigidity. The calculation is performed for a given amount 
of medium range order characterized by the fraction 77 of 
edge-sharing tetrahedra in the base glass {x = Q). This 
furtermore suggests that the strain can be concentrated 
in small rings and structures. 
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FIG. 2: Location of the rigidity (label 1) and the stress transi- 
tion (label 2) as a function of the equilibrium constant Ke for 
three fractions of edge-sharing tetrahedra in the base glass 
ri = 0.04 (solid line), rj = 0.33 (broken line) and rj = 0.50 
(dots). Vertical arrows serve to identify the width Ax of the 
intermediate phase at a given Ks- The insert shows the cor- 
responding width of the intermediate phase. 



\ X »• 

- \ c(2) 

~ \ stressed rigid 




1 

^ ^c(l) 


1 






--. K=2.5 

e 

K=0.1 

e 


- 












/ isostatic 











0.1 0.2 0.3 

Modifier concentration x 



FIG. 3: Probability of finding stressed rigid and isostatically 
rigid clusters with respect to the modifier concentration x for 
two different equilibrium constants: =0.1 (solid line) and 
Ks ~ 2.5 (broken line). Note that the stressed rigid cluster 
probability extrapolates to Xc ~ 0.20 in the MF description 
(bold line). The dotted vertical lines serve to define the two 
transitions and the intermediate phase for = 0.1. Here 
77 = 0.33. 



Fig. El shows the results of the construction. Large 
equihbrium constant (corresponding to a Q^-rich 
glass) will induce a large width for the IP. However, one 
observes that mostly affects the location Xc(2) of the 
stress transition whereas the location of the rigidity tran- 
sition Xf.(i) remains almost constant, as already signaled 
[lo| for a network glass. Noteworthy is also the shift from 
the MF rigidity transition at Xc = 0.20 to lower x that 
arises from the presence of weakly stressed rigid edge- 
sharing tetrahedra (EST, ric = 3.33 per atom). With a 
weaker rigidity due to the presence of these EST, part of 
the strain is captured in the EST, and onset of flexibil- 
ity can happen at lower modifier concentration x. As for 
IV- VI network glasses ilO|, the width Ax ~ Xc{2) — Xc{i) 
of the IP increases with the fraction rj of EST due to 
the shift of the location of the stress transition (Fig. [2J|. 
The change in speciation from a Q^-rich to a Q^-rich 
glass contributes however to an additional broadening of 
the IP. The trend with Ke observed in Fig. |2] can be 
further characterized from the computation of the prob- 
ability of stressed rigid and isostatically rigid clusters us- 
ing SICA. When the chemical reaction is desequili- 
brated on the left side (low Ke, i.e. a Q'^-rich glass), 
each modifier molecule will create mostly two flexible 
units (ric = 2.55 per atom) that serve to accumulate iso- 
statically rigid subregions of the network, as Q"* — 
connections are likely to appear. These are maximum at 
the stress transition (solid line. Fig. |3Jl, consistently with 
numerical simulations Q. On the other hand, a higher 



value of Ke leads to the growth of even more flexible units 
(Q^'s, He = 2.0 per atom) to the expense of Q'^'s, and 
will produce a stress transition at lower x (broken line. 
Fig. O. Indeed, increasing Ke at fixed x decreases the 
network mean coordination number and favours flexible 
— instead of isostatically rigid Q* — bondings. 
As a result, with growing concentration x, the network 
will lose stress earlier and will display a lower isostaticity 
in the IP. Thus Xe(2) is shifted to lower x. 

The constraint-related free eneregy is now considered, 
following the approach initially reported by Naumis .2^ . 
The free energy of the system is given by: 

Ke) = -/(^) + fesT^pg) Inpg^ (4) 

where — Z^'^-' is the stress energy equal to the number 
of redundant constraints, i.e. additional constraints that 
cannot be balanced by the degrees of freedom, and which 
vanish for x > Xe(i)- Figure 0] shows that the stress 
transition a,t x = Xe{2) is first order for any Ke- How- 
ever, with respect to the mean-field case (bold solid line) 
where Xe(i)) — Xc(2), the jump of the first derivative 
dJ-(2){x, Ke) / dx dX X = Xe{2) decrcascs with growing Ke- 
In the MF case, this jump is equal to 75.22, whereas 
it is only 45.30 for Ke ~ 2.5. This suggests that the 
transition broadens when the equilibrium |^ displaces 
to the right side, leading to a Q^-rich glass. On the other 
hand, the change in character with Ke of the rigidity 
transition is weak and second order. Finally, some chem- 
ical self-organization of the network is allowed through 
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mon. LPTMC is Unite Mixte de Recherche du Centre 
National de la Recherche Scientifique (CNRS) n. 7600. 
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FIG. 4: Free energy jr(2) (a;, ATe) — ^^(2) (0, ATe) of the system 
as a function of the modifier concentration x in the mean- 
field case (bold solid line), and for two different equilibrium 
constants: Ke — 0.1 (solid, line) and Ke ~ 2.5 (broken line), 
both for ri = 0.33. The arrows indicate the stress and rigidity 
transition for Ke =2.5. The insert shows the same quantity 
with X rescaled to its stress transition composition 2:c(2), to- 
gether with the free energy (dots) minimized by Ke- In all, 
ksT^ 1. 



an adaptative speciation. As stress costs energy, is is 
natural to imagine that the glass network will try to self- 
organize in the stressed rigid phase to decrease the energy 
by rewiring and reset some weaker bonds such as the ionic 
M-X ones that will lead to a Q" specie recombination. In 
the present description, this means that at a given con- 
centration X < Xc{2) 1 the minimization of the free energy 
can be accomplished with respect to Kf,. The equilib- 
rium constant minimizing J-{x, K^) provides then an 
estimation of the speciation via equ. (jSJl and the prob- 
ability of stressed rigid clusters. This leads to a stress 
transition at a; = Xc{2) — 0.168 and a free energy (dots 
in the insert of Fig. that is very close to the = 0.1 
speciation model and to sodium silicates 0| . 

In summary, with decreasing modifier concentration 
X there are different ways for a flexible system to self- 
organize in order to avoid stress, either by nucleating 
weak stress in small rings (EST), or by producing more 
flexible local structures (Q^'s) that balance the addition 
of new constraints arising from the decrease of the mod- 
ifier content. Both delay the onset of stressed rigidity. 
As a conclusion, we provide a prediction of the IP for 
sodium seleniosilicates (M=Na, X=Se) that have an EST 
fraction of 77 = 0.50 in the base network former SiSe2 
ji?! and an equilibrium constant = 0.15 [? ]. Ac- 
cording to the present approach, one therefore expects 
a stress transition at Xc{2) = 0.113, a rigidity transition 
at Xc{i) = 0.182 and a width for the Intermediate Phase 
of about Ax = 0.069. Larger structural correlations will 
probably refine this picture and are under consideration. 
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